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Abstract 

The low-energy scattering of ortho positronium (Ps) by H, He, Ne, and Ar atoms has 
been investigated in the coupled-channel framework by using a recently proposed time- 
reversal-symmetric nonlocal electron-exchange model potential with a single parameter 
C. For H and He we use a three-Ps-state coupled-channel model and for Ar and Ne we 
use a static-exchange model. The sensitivity of the results is studied with respect to the 
parameter C. Present low-energy cross sections for He, Ne and Ar are in good agreement 
with experiment. 
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Low-energy collision of ortho positronium (Ps) atom with neutral gas atoms and molecules 
is of interest in both physics and chemistry. Recently, there have been precise measurements 
of low-energy ortho-Ps scattering by H2, N2, He, Ne, Ar, C4H10, and C5H12 0, Due 
to internal charge and mass symmetry, Ps atom yields zero elastic and even-parity transition 
potentials in the direct channel. Ps scattering is dominated mainly by exchange correlation 
at low energies. If basis states are included for both Ps and target in a coupled-channel 
formalism, the number of channels grow as A^^. This complicates the tractability of the Ps 
scattering process in a coupled channel scheme compared to the electron scattering. The 
dominance of the short-range interaction in Ps scattering causes serious trouble towards the 
convergence of any coupled-channel formalism with a truncated basis |^. The use of 22 
coupled Ps pseudostates for Ps-H system in the R-matrix approach has indicated convergence 
difficulties []5|, ^ in Ps-H binding and resonance energies. 

To find a solution to the nonconvergence problem, we have suggested a nonlocal electron- 
exchange model potential 0, § with a single parameter C and demonstrated its effectiveness 
by performing quantum coupled-channel calculations using the ab initio framework of close 
coupling approximation. Two versions of this potential were suggested: one is time-reversal- 
symmetric and the other nonsymmetric. The nonsymmetric potential has been applied for the 
study of Ps scattering by H 0, He @] and H2 |§| using a three-Ps-state coupled- channel model. 
For He and H2 these studies yielded total cross sections in good agreement with experiments 
[|I|, 0, 1^ in addition to producing the correct pick-off quenching rate for Ps-He at low energies 
p. Higher Ps-excitation and ionization cross sections were also calculated in these cases to 
reproduce the total cross sections at medium and high energies. Previous theoretical studies 
10, 11 1 on Ps-He produced results in strong disagreement |Q with experiment [|l], ^ 

In a subsequent application of this model potential to Ps scattering by H |Q, it was found 
that the symmetric form leads to by far superior result than the nonsymmetric form. The 
symmetric form was able to reproduce accurate variational results very precisely for singlet 
Ps-H binding and resonance energies; the nonsymmetric form failed to yield such results. 

In view of this we reinvestigate the problem of low-energy Ps scattering by He 0] using the 
symmetric exchange potential employing the three-Ps-state model used before. We study how 
the low-energy cross sections for Ps-H and Ps-He change with the variation of the parameter 
C of the potential. Then we also apply this exchange potential to the study of low-energy Ps 
scattering by Ne and Ar using a static-exchange model. The present calculation accounts for 
the measured low-energy cross sections |1|, |^, |^ satisfactorily for He, Ne, and Ar. 

The total wave function of the Ps-target system is expanded in terms of the Ps eigenstates 
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as 



(ri,...,rN;rN+i,x) = ^ 



N 

+ 5](-l)^-+^'^F,(pj)x.(tj)0o(ri, FN+i, Tn) 



:i) 



where antisymmetrization with respect to Ps- and target-electron coordinates has been made. 
Here pi = (x + ri)/2, ti = (x — ri), where ri, ^ = 1, N denote the target electron coordinates, 
and fn+i and x are the electron and positron coordinates of Ps; 0o and Xu denote the target 
and Ps wave functions, and F^, is the continuum orbital of Ps with respect to the target. 
Sn+ij is the total spin of the Ps electron (A^ + 1) and target electron (j) undergoing exchange 
and can have values of 1 or 0. The spin of the positron is conserved in this process and 
the exchange profile of the Ps-target system is analogous to the corresponding electron-target 
system. Projecting the resultant Schrodinger equation on the Ps eigenstates and averaging 
over spin states, the resulting momentum-space Lippmann-Schwinger scattering equation for a 
particular total electronic spin state 5* can, in general, be written as ITT], IT2 



/,^,(k',k) = 05,,(k',k) 

1 \ " f ^^ /r k^O/^v(k^^ k) 



where /^j, is the scattering amplitude, and Bi^^i, is the corresponding Born potential, u and u' 



denote initial and final Ps states, k,yi' = y2m{S — e")/Ti? is the on-shell relative momentum of 
Ps in the channel v" with e" the total binding energy of the intermediate Ps and target states, 
and £ the total energy of the Ps-target system and m the reduced mass. Here 

N 

S5v(k', k) = 5,^(k', k) + E(-l)^""^'^^5t(k', k); (3) 

where is the direct Born potential and B^^ is the model potential for exchange between the 
Ps electron (denoted by + 1) with the jth target electron. For Ps-H scattering both S = 5*2, i 
= 0, 1 will contribute |jl2|; the corresponding potentials and amplitudes are usually denoted by 
B"^,,^ and /^^, respectively 0. For Ps scattering from He, Ne, and Ar etc. there will be only 
one scattering equation (Q) corresponding to total electronic spin 5* = 1/2. For these targets 
with doubly occupied spatial orbitals, in the sum over j in Eq. only half of the occupied 
target electrons in each sub-shell will contribute when the target is frozen to its ground state 



T3| . Consequently, only Sn+ij = 1 (/ and B ) will contribute to target-elastic scattering 



11 1 for targets with doubly occupied spatial orbitals. 
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For Ps-H scattering, the differential cross section is given by da/dQ = + 3|/^p]/4. 

For target-elastic Ps scattering by He, Ne, and Ar, the differential cross section is given by 
da/dQ = In all Ps-scattering the direct potential, B^, is zero for elastic and all even- 

parity-state transitions of Ps. Thus the nonorthogonal exchange kernel alone dominates the 
solution of Eq. and this dominance is possibly responsible for convergence difficulties to 
conventional approaches based on Eq. (^. 

The present exchange model was derived using Slater-type orbital for the H-atom, so that a 
generalization from a H-target to a complex target represented by a Hatree-Fock wave function 
becomes straight-forward. For Ps scattering from a H orbital, the model exchange potential 
between the Ps-electron (r2) and the orbital electron (ri) was derived from the l/r^ term and 
is given by [g If: 

4(_1)W f 

Bf^'u'^u = ^jj^ J 0;;'(r2)exp(zQ.r2)0;,(r2)dr2 

X I x:,(t2)exp(zQ.t2/2)x.(t2)dt2, (4) 

where / and /' are angular momenta of the initial (xi/) and final {Xu') Ps states, 0^ and 0^/ are 
initial and final H states, and Q = ki — kf . Here ki and kf are initial and final Ps momenta, 
respectively. In Eq. the symmetric form of the averaged quantity < D > is [H, ^] 



7,2 I 1.2 r^,2 I ^,2 

<D>= ^ + 



(5) 



where a^//2 and are the binding energies of the final target and initial Ps states, respectively, 
and C is the only parameter of the potential. Normally, the parameter C is taken to be unity 
0, ^ which leads to reasonably good numerical results. However, it can be varied slightly from 
unity to obtain a precise fit of a low-energy scattering observable (experimental or variational), 
as have been done in some applications of model potentials [p^ p!5[] . A variation of C from 
unity leads to a variation of the binding energy parameters (a^, etc.) used as average values 
for square of momentum |^ in the expression for (D) of Eq. (^. This, in turn, tunes the 
strength of the exchange potential (Q) at low energies. At high energies this model potential 
is insensitive to this parametrization and leads to the well-known Born-Oppenheimer form of 



exchange |T6|. We have turned this flexibility to good advantage by obtaining precise agreement 
with low-energy results of Ps scattering by H, He, Ne, and Ar, as we shall see in the following. 

For a complex target the space part of the HF wave function [|l^ is given by \I'(ri, r2, r^,- 
...,r]\r) = ^[0i(ri)02(i"2)---0j(i"j)---0Af(rAr)], where A is the antisymmetrization operator. The 
position vectors of the electrons are ri,i = 1,2, ...,iV and 0j's have the form: 0j(rj) = 
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Z^K c^Ki0ftj(i"j)- The orbital 0Kj(rj) is a Slater-type orbital. Considering proper antisymmetriza- 
tion with respect to Ps and target electrons, the final model exchange potential obtained from 
d) is given by [§ 

4aKjaK'j(-l)'+'' 



X / C'i(rj)exp(zQ.rj)0^j(rj)drj 

X J x:.(t,)exp(zQ.t,/2)x.(t,)dt,. (6) 



with 



where a^j is the energy parameter corresponding to the orbital 0Kj(rj) |[17|| . 

We use exact wave functions for H and Ps, HF wave functions for He, Ne, and Ar After 
a partial-wave projection, the one-dimensional scattering equations are solved by the method 
of matrix inversion. 

First we study the effect of the variation of the parameter C of the exchange potential 
in the Ps-H system using a three-Ps-state model with Ps(ls,2s,2p) states 0. We start our 
discussion with the singlet S-wave resonance. In Fig. 1 we plot the S-wave phase shift at 
different energies which illustrate the resonance pattern obtained with different values of C. 
The resonance position shifts monotonically towards lower energies with decreasing of C from 
unity. We have shown it in steps where the value of C is varied from unity to 0.785. The 
resonance position matches with the accurate prediction of 4.01 eV [1^ for C = 0.785. In Fig. 
2, we plot kcot6 versus k"^ for the corresponding low-energy S-wave phase shifts 6. Figure 2 
demonstrates how the improvement in the resonance position simultaneously improves the Ps-H 
binding energy. For C = 1.0 the resonance and binding energies are 4.715 eV and 0.165 eV, 
respectively; for C = 0.9 the corresponding energies are 4.470 eV and 0.445 eV, respectively. 
At C = 0.785, while the resonance position is correctly fitted to 4.01 eV (Fig. 1), we obtain 
an approximate binding energy of 1.02 eV from a linear extrapolation as in Fig. 2, and 0.99 
eV with more precise fitting considering second order corrections, compared to the accurate 
prediction of 1.0598 — 1.067 eV [|l^. This behavior of the low-energy phase shifts, which yields 
simultaneously the Ps-H resonance and binding energies, indicates that the use of the model 
potential in a coupled-channel scheme can lead to a good description of Ps-H scattering. 

We exhibit in Figs. 3 (a), (b), and (c) the present elastic cross sections, for Ps scattering 
by He, (three-Ps-state [|]) Ne and Ar (static-exchange), respectively, for C = 1, 0.85, and 
0.785. The value C = 0.785 yielded the good agreement in the case of Ps-H. For the closed- 
shell atoms, a typical C close to 0.85 works well for the 3-Ps-state model in Ps-He and for 
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the static-exchange model in Ps-Ar and Ps-Ne. Although the present cross sections differ from 
other theoretical |T^, [11]] and experimental works at low energies for Ps-He (See, Fig. 6 of 
Ref. 1^), they agree well with the recent measurements of Skalsey et al. 0] and unpublished 
work of G. Peach as quoted in Ref. 0. 



Table I: Low-energy S-wave phase shifts in radians for Ps-He, Ps-Ne, and Ps-Ar for different 
k in au. The entries for k = correspond to the scattering lengths in units of oq, incident 
positronium energy E = G.Sk"^ eV. 



k 


Ar(SE) 


Ne(SE) 


He(SE) 


He(3st) 


0.0 


1.65 


1.41 


1.03 


0.90 


0.1 


-0.164 


-0.141 


-0.103 


-0.088 


0.2 


-0.319 


-0.277 


-0.202 


-0.172 


0.3 


-0.457 


-0.404 


-0.294 


-0.249 


0.4 


-0.572 


-0.518 


-0.375 


-0.315 


0.5 


-0.656 


-0.615 


-0.444 


-0.368 


0.6 


-0.706 


-0.694 


-0.500 


-0.408 


0.7 


-0.720 


-0.754 


-0.541 


-0.433 


0.8 


-0.699 


-0.792 


-0.569 


-0.445 



Table II: Low-energy P-wave phase shifts in radians for Ps-He, Ps-Ne, and Ps-Ar for different 
k in au. The incident positronium energy E = 6.8/c^ eV. The numbers in parenthesis denote 
powers of ten. 



k 


Ar(SE) 


Ne(SE) 


He(SE) 


He(3st) 


0.1 


-2.71(-3) 


-1.63(-3) 


-8.19(-4) 


-6.24(-4) 


0.2 


-1.94(-2) 


-1.20(-2) 


-6.11(-3) 


-4.63(-3) 


0.3 


-5.51(-2) 


-3.55(-2) 


-1.85(-2) 


-1.38(-2) 


0.4 


-1.06(-1) 


-7.17(-2) 


-3.81(-2) 


-2.80(-2) 


0.5 


-1.65(-1) 


-1.17(-1) 


-6.33(-2) 


-4.53(-2) 


0.6 


-2.23(-l) 


-1.66(-1) 


-9.15(-2) 


-6.26(-2) 


0.7 


-2.75(-l) 


-2.15(-1) 


-1.20(-1) 


-7.63(-2) 


0.8 


-3.17(-1) 


-2.60(-l) 


-1.47(-1) 


-7.89(-2) 



In Table I we present the S-wave phase shifts for Ps-Ar (static exchange model denoted 
SE), Ps-Ne (SE), and Ps-He (SE and three-Ps-state models) scattering for C = 0.85. In Table 



6 



II we present the same for the P wave. The magnitude of the scattering lengths, low-energy 
cross sections and phase shifts (well below Ps-excitation threshold) increase monotonically as 
we move from He to Ne and from Ne to Ar. As the effective potential for elastic scattering 
in these cases is repulsive in nature, this indicates an increase in repulsion from He to Ne and 
from Ne to Ar. 

We find from Figs. 3 that the energy-dependences of the elastic cross sections are similar 
for all the closed-shell atoms studied here. The cross section has a monotonic slow decrease 
with increasing energy. This trend is consistently found in all previous theoretical calculations 
in Ps-He. Also, this is expected as the underlying effective potential for elastic scattering is 
repulsive in nature. 

In conclusion, we have reinvestigated the problem of low-energy elastic Ps scattering by 
H, and He (three-Ps-state) using a symmetric nonlocal electron-exchange potential with a 
parameter C. We further apply this potential to Ps scattering by Ne, and Ar (static-exchange). 
Although the value C = 1 was originally suggested, a slightly lower value C ~ 0.8 leads to 



good agreement with accurate experiment and accurate calculations ||T8|, |I9| in the present 
cases. Although, a non-symmetric form of the model potential provides a fairly good account 
of the cross section [0 , we have found that the symmetric form is able to provide a more precise 
description of scattering. The Ref. we have demonstrated the effectiveness of the present 
exchange potential in electron-impact scattering. Simplicity of the present exchange potential 
and the reliability of the present results calculated with it from a two- (Ps-H) to a 19-electron 
(Ps-Ar) system reveals the effectiveness of the exchange model and warrants further study with 
it. 

We thank Prof. B. H. Bransden for his helpful and encouraging comments. The work is 
supported in part by the CNPq and FAPESP of Brazil. 
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Figure Caption 

1. S-wave singlet Ps-H phase shifts in radian showing the variation of the resonance position 
with the variation of C in (|^) using present three-Ps-state model. 

2. k cot 6 and ik versus k'^ plot showing the change in Ps-H binding energy with the variation 
in C as in figure 1 (Energy = 6.8A;^ eV). The crossing of the k cot 6 and ik curves give the energy 
of the bound state. 

3. Cross section of Ps scattering by (a) He, (three-Ps-state model) (b) Ne, and (c) Ar 
(static-exchange model) for different C: C = 1 (full line), C = 0.85 (dashed-dotted) line, 
C = 0.785 (dashed line), experiment (box, Ref. @). 
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Figure 1 
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Figure 2 
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